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Abstract. The nature of dark matter is still unknown and one of the most fundamental 
scientific mysteries. Although successfully describing large scales, the standard cold dark 
matter model (CDM) exhibits possible shortcomings on galactic and sub-galactic scales. It 
is exactly at these highly non-linear scales where strong astrophysical constraints can be set 
on the nature of the dark matter particle. While observations of the Lyman-a forest probe 
the matter power spectrum in the mildly non-linear regime, satellite galaxies of the Milky 
Way provide an excellent laboratory as a test of the underlying cosmology on much smaller 
scales. Here we present results from a set of high resolution simulations of a Milky Way 
sized dark matter halo in eight distinct cosmologies: CDM, warm dark matter (WDM) with 
a particle mass of 2 keV and six different cold plus warm dark matter (C+WDM) models, 
varying the fraction, /wdm? and the mass, mwdm? of the warm component. We used three 
different observational tests based on Milky Way satellite observations: the total satellite 
abundance, their radial distribution and their mass profile. We show that the requirement of 
simultaneously satisfying all three constraints sets very strong limits on the nature of dark 
matter. This shows the power of a multi-dimensional small scale approach in ruling out 
models which would be still allowed by large scale observations. 
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1 Introduction 

On large scales cold dark matter (CDM) predictions agree very well with galaxy clustering 
surveys such as the Sloan Digital Sky Survey (SDSS) and 6dF. On galactic and sub-galactic 
scales some tension between observations and numerical CDM simulations appears and hence 
there are a some small scale challenges for the CDM paradigm. In the following we want to 
emphasize three issues discussed in the literature. The first one is the so called core/cusp 
problem: the estimated mass profiles of dark matter haloes inferred from observations are 
mostly consistent with constant density cores [e.g. 1,2], whereas simulations of CDM haloes 
tend to produce cuspy inner density profiles [e.g. 3-5]. The second concern is the overabun- 
dance of substructure in Milky Way sized dark matter haloes when compared with observed 
satellite galaxies [6, 7]. A third issue is that several of the most massive CDM subhalos have 
higher central densities than the observed in any of the satellite galaxies. The ACDM model 
predicts the existence of large, dark galactic subhalos, which seem to be "too big to fail" to 
form a dwarf galaxy [e.g. 8-11]. 

Even though several baryonic solutions have been proposed to solve at least some of the 
inconsistencies mentioned above [e.g. 12], it is worth exploring alternative scenarios, where 
the dark matter density of the Universe is not made of a cold and thermal relic. If the dark 
matter particle is much lighter (with a mass in the keV region) than its CDM counterpart, 
it decouples from the hot plasma in the early Universe when it is still relativistic. Hence, 
dark matter particles are able to free-stream until relatively late times, washing out small 
density perturbations and suppress the formation of structure below a characteristic scale 
(which is denoted as the free streaming scale, Afg). Such dark matter models are named 
Warm Dark Matter (WDM) [e.g. 13, 14, and references therein]. Typical candidates are 
the supersymmetric gravitino [15-17] and the non-thermally produced sterile neutrino [e.g. 
18-20]. 

Another viable solution is to assume that the total dark matter density is made of a 
mixture of cold plus warm dark matter (C+WDM) [21-24]. In addition to the mass of the 
warm particle, these class of models are also characterised by the fraction of warm to cold 
dark matter, i.e. 

_ ^WDM _ ^WDM ^-^ -^^ 

^WDM + ^CDM ^DM 
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Although afflicted with large uncertainties, currently, the best way to distinguish between 
CDM and non-standard models is to analyse data from the Lyman-a forest [for a review see 
25]. As we will show in this work, a dark matter model, even when not ruled out by Lj-a 
physics, still has to pass a series of test on much smaller scales: in order to be hosts of the 
brightest Milky Way satellites, subhalos in a galaxy-sized dark matter halo should be in line 
with observations in terms of their abundance, radial distribution and mass profiles. 

This paper is organised as follows. Section 2 provides a brief summary of the techniques 
we used for computing the transfer function for mixed dark matter models as well as an 
overview of our numerical simulations (see also [26] for more details). Section 3 shows the 
results of our subhalo analysis: subhalo abundance, radial distribution and mass profiles are 
compared with observed satellite galaxies. A discussion and conclusion is presented in section 
4. 

2 Cosmological Models and N-body Simulations 

In this section we briefly introduce the cosmological models used in this work, as well as 
the simulation details. Besides the conventional cold dark matter paradigm, we run a warm 
dark matter cosmology with a particle mass of 2 keV as well as six different cold plus warm 
dark matter (C+WDM) cosmologies. It should be emphasized that all non-standard models 
except two (f20 - O.lkeV & f50 - 0.3keV) lie within the two sigma range of validity for Lj-a 
forest constraints [see e.g. 22, for recent constraints]. 

A generic C+WDM model is described by the mass of the warm particle, mwDM, and 
the fraction of warm to cold dark matter, /w (see eq. (1.1)). We have chosen a set of six 
mixed dark matter models by scanning through the /w - ^wdm parameter space. In order 
to isolate effects coming from either the fraction or the mass, in some of the models only one 
of the parameters is changed. Unlike the numerical C+WDM realisations in [22, 24], where 
every particle in the simulation represents the mixture of both species, our two particle 
species approach allows us to distinguish the cold and warm components (similar as for 
example in SPH simulations containing baryons). The initial conditions are set up using two 
cosmological boxes, one containing only cold particles, the other only warm particles, and 
merging them with an offset of half a grid cell in each dimension. A detailed description 
of the numerical methods is provided in section 2 of [26], hence in the following we just 
briefly sketch the simulation procedure and the implementation of the corresponding transfer 
functions. At first, we run large scale simulation boxes with length of 40 Mpc with a total 
particle number of 2 x 256^. The transfer functions of the non-CDM models, representing 
small scale suppression, can be written as 

T\k) = 7^^^, a = {WDM, C+WDM}. (2.1) 

For the pure WDM run, we used the fitting formula suggested by [27], 

Twdm(A:)' = (l + (aA:)'")"''^^ (2.2) 

where v — 1.12 and a is a function of the warm dark matter particle mass. In the mixed 
dark matter case, we made the simplistic assumption that the universe is filled with a dis- 
sipationless, coupled cold plus warm dark matter fluid and numerically solved the basic 
hydrodynamical equations of Newtonian physics. Linear theory then allows us to connect 
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Figure 1. Left panel: set of N-body simulations in the /w - ^wdm parameter space. Green stars 
denote mixed dark matter simulations whereas red stars are pure cold (/w = 0) and pure warm 
{fw = 1) dark matter simulations respectively. Right panel: matter power spectrum of all simulated 
models with increasing /w from top to bottom (at large k). 



Table 1. Details of the simulations. i^200, ^200 and M200 are measured with respect to 200 times 
the mean matter density. 



Label 


mwDM 


fw 


z=99 


-R200 


-^200 


M200 




[keV] 


[%] 


[km s"-*-] 


[kpc] 


[xlO^] 


[IO12 Mo] 


CDM 










367 


9.88 


1.36 


f05 - 0.05keV 


0.05 


5 


82.1 


362 


8.98 


1.24 


f05 - O.lkeV 


0.1 


5 


32.6 


368 


9.46 


1.31 


f20 - O.lkeV 


0.1 


20 


51.7 


365 


8.49 


1.17 


£20 - O.SkeV 


0.3 


20 


11.9 


368 


9.20 


1.27 


£50 - O.SkeV 


0.3 


50 


16.2 


367 


9.01 


1.24 


£80 - IkeV 


1.0 


80 


3.8 


368 


8.69 


1.20 


WDM 


2.0 


100 


1.6 


367 


9.04 


1.25 



density contrasts between two epochs (i.e. matter-radiation equality and the starting redshift 
of the simulations) via a mixed dark matter transfer function T{k) [see 26, for details]. The 
right panel in figure 1 shows the linear power spectrum at the starting redshift of our simula- 
tions zic = 99. Unlike the pure WDM case, where the linear power rapidly drops to zero, the 
linear power spectrum in mixed dark matter models stabilizes due to the presence of a cold 
component and approaches a constant plateau with a characteristic height only dependent on 
the fraction /w- All simulations have been performed with the parallel treecode PKDGRAV, 
written by Joachim Stadel and Thomas Quinn [28]. Initial conditions are generated with a 
parallel version of the GRAFIC package [29] and are based on the cosmological parameters 
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Figure 2. Projected density maps of the refined Mill^y Way sized halo, varying the cosmology. Each 
box shows a sphere with radius r2oo- Top row from left to right: CDM, WDM 2keV, f05-0.1keV, 
f20-0.3keV. Bottom row from left to right: f80-lkeV, f05-0.05keV, f50-0.3keV, f20-0.1keV. 



taken from WMAP7 [30]: ag = 0.8, h = 0.7, f^dm = 0.227, f^b = 0.046, ft a = 0.727 and 
Us = 0.961. For the sake of completeness, figure 1 shows the resulting grid of high resolution 
N-body simulations. 

In a second step we selected an isolated and Milky Way sized IO^^Mq) halo at 
redshift z = and re-simulated a Lagrangian region of 3 x i?200 of the original object at 8^ 
times better mass resolution. This leads to a particle mass of rrip = 1.38 x IO^Mq and a 
gravitational softening of 355 pc in the refined region. 

Due to a decoupling when still relativistic, WDM particles are expected to retain a 
thermal motion that is assumed to influence the internal structure of dark matter haloes [e.g 
5, 31, and references therein]. We therefore added a thermal velocity component to all WDM 
particles on top of the Zeldovich velocities at the starting redshift zic = 99, by following the 
method presented in refs. [5, 26]. An overview of all six mixed dark matter simulations and 
the corresponding pure CDM and pure WDM runs can be found in table 1; the projected 
density maps are shown in figure 2. 

Finally, subhalos in all refinement runs are identified using 6dFoF [32], a phase-space 
extension of the well known FoF method. The minimum number of particles in a subgroup 
is set to A^min = 64. The subhalo properties such as mass and maximum circular velocity are 
computed with a profiling routine that takes the phase-space coordinates of 6dFoF as input. 
Details can be found in [33]. 

3 Subhalos in C+WDM 

Dark matter subhalos have implications for both dark matter detection and satellite galaxy 
properties [see e.g. 34, 35, for recent reviews]. Here we present the abundance, spatial distri- 
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bution and mass profiles of subhalos. Comparing all of these quantities with observations of 
satellite galaxies provides a solid test of the underlying cosmology. 

Some models of galaxy formation connect the most massive subhalos at redshift z — 
with the most luminous dwarf galaxies observed in the Milky Way [e.g. 36, 37]. More 
recent studies suggest that the observed Milky Way satellites correspond to the most massive 
subhalos at the time of accretion [e.g. 38-40]. Caused by the gravitational potential of the 
main halo, subhalos experience tidal stripping, which can lead to a substantial mass loss of 
the subhalo after infall. 

Instead of halo and subhalos masses we use the maximum circular velocity, i^max^? which 
is a well defined proxy for halo and subhalo size and easy to determine accurately from cos- 
mological N-body simulations (both at z = and at infall). Extracting the maximum circular 
velocity at the time of accretion (i^peak = '^max( ^infall)) foi" each halo is computationally more 
expensive than just taking the corresponding value at redshift z = (i.e. i;max(^ = 0)), since 
it requires the assembly of a merger tree. 

In appendix A we present a simple method how to compute an estimate for i;peak — 
'^max( ^infall) of ^ subhalo, bascd on its present-day properties alone (i;max(^ = 0) and binding 
energy). Although afflicted with some scatter, it turns out to be a useful estimator for 
the maximum circular velocity at infall and provides quick method for a comparison with 
observational data. For the following analysis, we will make use of two selection criteria, 
'^max and i^peak? whenever comparing with observational data. 

3.1 Abundance 

In the past years, the SDSS survey has lead to the discovery of more than 10 very faint 
Milky Way satellites [e.g. 41, 42]. This alleviates the original missing satellite problem [6, 7] 
somewhat, but there is still a difference up to half an order of magnitude between simulations 
and observations [e.g. 10, 46-48]. The left panel of figure 3 shows the subhalo mass function, 
the right panel shows the cumulative velocity function for all subhalos with a value of i;max > 9 
km s~^ (which is approximately the resolution limit of the CDM reference simulation), as 
well as a sample of observed Milky Way satellites (black open circles). When representing 
the subhalo abundance as a function of tidal mass, most of the non-standard models closely 
follow the CDM case. On the other hand, plotting the subhalo abundance as a function of 
'^max5 the differences get larger since i;max is sensitive to the concentration of the subhalo. 
The observation sample consists of both classical and SDSS satellites, though for the latter, 
we included a completeness correction since the SDSS only covers about 28 per cent of the 
sky. Whereas the CDM cumulative velocity function is well approximated by a power law 
[e.g. 49, 50], 

'^max) ^ ('^max/ '^max,host (3.1) 

all cosmologies containing some fraction of warm dark matter exhibit a reduced total abun- 
dance of substructure (and therefore a shallower velocity function) due to the lack of power 
on small scales. As shown in [26], a generic mixed dark matter cosmology containing a ther- 
mal WDM candidate can be brought into a hierarchy according to its effective free-streaming 
length, i.e. 

^The maximum circular velocity is defined at that radius rmax at which — y^GM{< r)/r gets maximised. 
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Figure 3. Subhalo abundances in all cosmologies, represented as a function of subhalo tidal mass 
(left panel) and as a function of the maximum circular velocity Vmax at redshift z = (right panel). 
Estimates for the Vmax of the host subhalos of the Milky Way satellite galaxies are shown for compar- 
ison (open circles, right panel) [8, 43, 44]. Note that converting from the observed stellar kinematics 
to the much larger scale of I'max introduces large uncertainties and relies on CDM subhalo scaling 
relations [45]. 



Hence, the C+WDM models f20 - O.lkeV and f50 - O.SkeV are the most extreme in our 
sample. This is not only apparent in the linear power spectrum (right panel of figure 1) but 
also in the subhalo velocity functions and mass functions (orange and green solid lines): only 
a few subhalos exist in galactic halos for these two models. 

Figure 3 also shows the estimated host subhalo sizes for the observed Milky Way satellite 
galaxies. Since the host subhalo rmax lies well beyond the galaxies, these size estimates are 
very uncertain and it is more accurate to compare simulated subhalo mass profiles with 
observations at much smaller radii [45] , especially the half light radii [44] , which we will do 
in section 3.3. Using the approximate satellite galaxy i;max values for now, we find that pure 
2keV WDM scenario as well as the mixed dark matter models f80 - IkeV, f20 - O.SkeV and 
f05 - 0.05keV reproduce the observed abundance quite well. In the pure CDM and f05 - 
O.lkeV case, the classical problem of missing satellites is clearly visible and some mechanism 
to suppress star formation in some fraction of these subhalos is required to match the velocity 
function estimated from the observations. It is possible that with increasing sensitivity of 
observational methods over the next several years, many more faint and ultra-faints dwarfs 
will be detected, which might continuously close the gap between the number of observed 
dwarf galaxies and simulated CDM subhalos [e.g. 51]. For a very solid, conservative constraint 
on the underlying cosmology we therefore just demand that the number of currently know 
satellites be comparable or larger than the number of subhalos in the corresponding size range. 
This lower bound on the required subhalo abundance clearly rules out the two extreme mixed 
models f20 - O.lkeV and f50 - 0.3keV. 
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Figure 4. Top: Normalised radial distribution of all subhaloes within R200 that have that have 
'^max > 9 km/s (left) and 'L'peak >13 km/s (right) respectively. Bottom: Normalised radial distribution 
of the 12 subhalos with the largest values of the maximum circular velocity, selected after ^'max (left) 
and '^peak (right). Moreover, the normalised mass distribution of the main CDM halo is shown in 
both panels (black dashed line). The open circles refer the 12 most luminous satellites observed. 
Horizontal error bars refer to observational uncertainties. Starting with the smallest distance those 
are: Sagittarius, LMC, SMC, Ursa Minor, Draco, Sculptor, Sextans, Carina, Fornax, Leo II, Canes 
Venatici I, Leo I (data taken from [8, 43, 44]). 



3.2 Radial Distribution 

The radial distribution of satellite galaxies provides another important observational con- 
straint: it needs to be matched by any realistic set of subhalos selected from a cosmological 
simulation in order to be possible hosts of the satellite galaxies. The free streaming effect in 
warm dark matter models does not only suppress the subhalo abundance, but they are also 
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found to have later formation times and therefore lower concentrations when compared with 
CDM [e.g. 9, 47]. Lower concentration subhalos are more affected by tidal mass loss and less 
likely to survive in the inner halo, where some of the satellite galaxies are observed. 

Figure 4 shows the normalised subhalo radial distribution of all six mixed dark matter 
models as well as the two pure CDM and WDM runs. In the two panels on top, all subhalos 
with maximum circular velocity Vma^ > 9 km/s (left) and i^peak > 13 km/s (right) are 
used. The bottom panels only shows the 12 largest objects according to each of the two 
selection criteria used above (i;max on. the left and i^peak on. the right) as well as the 12 most 
luminous satellites observed. These consist of 11 classical dSphs (Sagittarius, LMC, SMC, 
Ursa Minor, Draco, Sculptor, Sextans, Carina, Fornax, Leo II, Leo I) plus Canes Venatici I 
with a luminosity comparable with Draco but only observed in the SDSS survey. The black 
dashed line refers to the normalised mass distribution of the CDM halo. 

As already mentioned, there is a theoretical motivation to believe that 'I'peak correlates 
better with galaxy luminosity than 'Umax- We will therefore mainly focus on the bottom 
right panel of figure 4. For small galactocentric distances (r < 100 kpc), the two models 
containing only a small WDM fraction as well as the pure 2 keV WDM model match the 
observed satellite distribution remarkably well. On the other hand, in mixed cosmologies 
with larger free-streaming length, i.e. f20 - O.lkeV, f20 - O.SkeV and f50 - O.SkeV, none of 
the 12 most massive subhalos reside on orbits close (r < 100 kpc) to the galactic center. 

3.3 Mass Profiles 

It was recently shown that internal structure and abundance of the bright Milky Way satellites 
is not consistent with populating the largest CDM subhalos in main halos of Myir ^ 1.7— 2.0 x 
10^^ M0 [8, 10, 52]. To match the observations some of the largest, densest CDM subhalos 
would have to fail to form a luminous dwarf galaxy, which seems rather unlikely. This issue 
has been quite accurately and very timely referred to as the "too big to fail" problem [8]. 
However, it turns out that most but not all of the halos in this mass range do contain too 
many large subhalos [e.g. 53], which raises the question why the Milky Way halo should be 
special. Since the discrepancy in the Via Lactea II simulation is only about a factor of two in 
mass [10] rescaling this halo to a smaller virial mass of (M - 0.8 X 10^2 Mq) would resolve the 
discrepancy. However, recent studies based on kinematic measurements of observed satellites 
[e.g. 54-56], as well as the presence of two massive objects as the Magellanic clouds [57] 
disfavour such a low halo mass for the Milky Way. Other ways to reconcile CDM with 
these constraints are discussed in the literature: For example baryonic feedback processes or 
highly stochastic galaxy formation in the most dense subhalos but it is currently unclear how 
plausible they really are [e.g. 52, and references therein]. 

Non-standard cosmologies like self-interacting, warm or mixed dark matter might pro- 
vide a solution to this issue and they have recently attracted a lot of interest [e.g. 9, 26, 58]. 
Lovell et al. (2012) [9] have re-run the Aquarius A halo in a thermal 1.4 keV warm dark 
matter cosmology. Indeed, they find that subhalos are not only less numerous but also less 
concentrated (due to their later formation time) than their CDM counterparts. Therefore, a 
WDM cosmology with a thermal dark matter candidate with a mass of roughly 2 keV does 
not suffer from the "too big to fail" problem. 

Figure 5 shows the circular velocity profiles at redshift z = for the 12 subhalos with 
the largest value of i^peak ioi all cosmologies under consideration. As mentioned above, it is 
evident that CDM predicts several subhalos which are too massive in order to be possible 
hosts of the most luminous satellites. Due to the reduced power at small scales, structure 



-8- 




r [kpc] r [kpc] 

Figure 5. Circular velocity curves at redshift z = for the 12 subhaloes that had the highest value 
of '^peak- Measured '^0110(^1/2) values for the brightest dwarf galaxies of the Milky Way are plotted as 
data points with error bars. In contrast to Fig. 4, LMC, SMC and Sagittarius are removed from the 
sample. 

formation in all other cosmologies is delayed, leading to slightly smaller halo and subhalo 
concentrations. When going towards models with increasing effective free streaming length 
(from left to right and from top to bottom), the central masses of the subhalos are substan- 
tially lowered. While the pure 2 keV WDM, f05 - 0.05keV and f80 - IkeV models match the 
observed satellite galaxies quite well, the remaining cosmologies, especially the two outlying 
models f50 - O.SkeV and f20 - O.lkeV, fall into the other extreme: they do not produce enough 
massive subhalos in order to host the observed satellite galaxies. 
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4 Discussion & Conclusions 

So far, the best observational constraints on the non-hnear matter power spectrum, and 
therefore on the nature of dark matter are given by the Lyman alpha forest (Ly-a) [e.g. 27]. 
Assuming a pure warm dark matter cosmology, Ly-a analysis provide a lower bound on the 
dark matter particle of roughly the order of one keV (strongly dependent on the production 
mechanism). In this paper we have demonstrated how satellite galaxies of the Milky Way 
with masses IC^^"^^ Mq allow to test dark matter on even smaller scales than the Ly-a 
constraints, i.e. we can probe (and rule out) free-streaming scales, which lie significantly 
below the Ly-a constraints. Therefore, a detailed comparison between observations and 
subhalos in a numerical simulation is a further powerful test to constrain the nature of dark 
matter. 

We have performed a set of of high resolution N-body simulations of a Milky Way sized 
dark matter halo in a variety of cosmological models: CDM, WDM with a particle mass of 2 
keV as well as six different mixed dark matter models (C+WDM) varying the fraction and 
the mass of the warm component (see left panel of Fig. 1 for an overview). The most massive 
subhalos in each simulation had then to pass three tests in order to be considered as a viable 
cosmological model. They should simultaneously reproduce the abundance of satellites, their 
radial distribution and their mass profile. Our findings can be summarised as follows: 

• Comparing the subhalo abundance with observations only provides a lower bound on a 
given cosmology in terms of its free streaming length. The number of subhalos in the 
models f20 - O.lkeV and f50 - 0.3keV hes clearly below the observed satellite galaxy 
abundance and therefore these models are ruled out. CDM and f05 - O.lkeV offer 
more than enough subhalos and could match the observational counts if some strong 
feedback mechanism suppresses galaxy formation in small halos, or in other words, 
both these models may suffer from the well known "missing satellites problem". The 
three remaining C+WDM models (f05 - 0.05keV, f20 - 0.3keV & f80 - IkeV) as weh as 
the pure 2keV WDM model are in good agreement in the range i^max ^ 20 km/s, i.e. 
they would produce a reasonable number of satellite galaxies, if each subhalo contains 
a galaxy. 

• Considering also the radial distribution of subhalos, two additional C+WDM cosmolo- 
gies are at the borderline of being ruled out. The largest subhalos (selected at infall 
time) of the models f20 - O.SkeV and f80 - IkeV reside almost entirely in the outer 
regions of the main halo, while the satellite galaxies tend to follow the main halo mass 
distribution. Also note that subhalo samples selected by their present (z=0) size are 
ruled out in all cosmologies (including CDM), because such samples are far too extended 
to match the observations. 

• Incorporating also the mass profiles into the big picture, it seems that only the pure 
WDM model with a 2keV particle is able to match the all observations. The CDM 
halo and the f05 - O.lkeV halo contain several subhalos that are clearly too dense in 
order to host the brightest satellite galaxies, i.e. they both suffer from the 'too big to 
fail' problem, but for the f05 - O.lkeV the discrepancy is much smaller. Model f80 - 
IkeV on the other hand, does quite well in the mass profile test, but it predicts a too 
extended spatial distribution of satellite galaxies. 

We can conclude that our detailed comparison between observations and subhalos rules out 
a range of mixed C+WDM models (f20 - 0.3keV, fSO - IkeV and f05 - 0.05keV), which are 
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not excluded by the Lj-a constraints. The most stringent test comes from the subhalo mass 
profiles (figure 5), which rules out four models because their subhalos are not dense enough; 
and it challenges CDM, which has several subhalos which are too dense. One model (f80 
- IkeV) gives an acceptable match in this test, but fails to match the radial distribution 
of observed satellites. This illustrates that for the best, strongest constraints one has to 
take into account both the densities and the radial distributions. The two tests are not 
independent however: most models, which do not form dense enough subhalos (f20 - 0.3keV, 
f05 - 0.05keV, f50 - 0.3keV and f20 - O.lkeV), also do not predict enough satellite galaxies in 
the inner halo. That is not surprising, since high central densities are needed for subhalos 
to resist tidal stripping in the inner halo [e.g. 59]. Furthermore, it is fair to mention that 
halo-to-halo scatter could also somewhat relax our conclusions. 

Whether or not baryonic physics will alter these conclusions is still an open question, 
since feedback processes might be able to reduce the dark matter densities near the centres of 
halos and subhalos. However, dwarf galaxies are known to be the most dark matter dominated 
objects in the Universe and it is therefore reasonable to assume that spatial distribution and 
mass profiles of these objects are mainly driven by the dark component and not by baryons. 
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A Correlation between orbital energy and subhalo mass loss 

As pointed out in the Introduction, most recent models of galaxy formation connect the most 
luminous satellites with the largest subhaloes at infall, i.e. the relevant subhalo size is the one 
measured at their first infall into the main halo (^^infaii) and not the present {z = 0) size that is 
left over after tidal mass loss, which can be very large for some subhalos (especially for those 
near the host halo center). The size at infall can be quantified with i;peak — '^max (^infall) ^- This 
requires for each subhalo in a simulation a merger tree built from a large number of simulation 
outputs. It has recently been shown that there is a strong correlation between the infall time 
of a subhalo and its binding energy at redshift z = [60, 61]. In the following, we present a 
strong correlation between the z=0 subhalo orbital energy and the ratio i;max(^ = 0)/i;peak- 
This provides a new and simple method how the required information, i.e. Vpeak^ can be 
extracted from only present day subhalo data, without building a computationally expensive 
merger tree. 

For this subhalo analysis, we use the Via Lactea II (VL2) simulation^ [49], a high 
resolution ACDM realisation of a refined dark matter halo with a mass similar to the Milky 

^For simplicity, ^;max(2;infaii) is denoted by ^^peak- This notation is motivated by the observation that most 
subhalos reach the peak vaule of their maximum circular velocity at infall [59] . 

^The subhalo data set of VL2 simulation can be found here: www.itp.uzh.ch/ diemand/vl/data.html. 
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Figure 6. Relation between the orbital energy of each subhalo of the VL2 sample and the ratio 
between Vma.^{z = 0) and 'ymax(^infaii), denoted by Vpeak- 



Way. We define the binding energy of a subhalo as the sum of kinetic plus potential energy, 
divided by the mass of the subhalo [62], 

eorbit = K -\-U 

where v is the velocity of the subhalo with respect to the halo center and the main halo is 
fitted with an NFW mass distribution [63] with scale radius of 21.5 kpc. A subhalo is 
considered to be bound to the main halo if eorbit < 0. Furthermore, only subhalos are listed 
in the sample that (a) lie within the virial radius of the main halo and (b) have a maximum 
circular velocity at infall of 'ymax(^infaii) ^ lOkm/s. This reduces the original set to a sample 
of 953 objects. 

Fig. 6 shows the median values of the ratio i;max(^ = 0)/'Upeak ioi these 953 subhalos 
binned by their orbital energy (red solid dots). The ratio Vma^iz = 0) /'Umax (^infall) is a proxy 
for the subhalo mass loss since the time of accretion due to tidal stripping. The green shaded 
area shows the region between the 16*^ and 84*^ percentiles. A tight correlation is clearly 
visible: the deeper a subhalo orbits within the potential well of the main halo, the more tidal 
stripping it has undergone. This leads to a reduced value of the maximum circular velocity 
at present time relative to the peak value at infall. The black solid line in Fig. 6 is a simple 
regression with respect to the median, given by the following one parameter fitting form: 

'^max 



'^peak 



(eorbit) = 1.0 + ^ • eorbit- (A.2) 
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Figure 7. Left panels: circular velocity curves fo the 12 most massive subhalos at redshift z = 
('^max(^ = 0))- Right panels: circular velocity curves of the 12 most massive subhalos at infall ('^peak), 
calculated with our correction scheme (Eq. (A. 2)). Each is illustrated for the pure CDM (top) as well 
as the pure WDM (bottom) simulation. The data points refer to the 9 brightest satellite galaxies in 
the Milky Way. 

The value of the fitting parameter is obtained using the Levenberg & Marquart method, 
yielding ^ = 2.89 • 10~^ [km/s]~^. The boundary condition of the above fitting form was 
chosen such that unbound subhalos (corbit = 0) have no mass loss {v 

max — '^peakj- Equation 

(A. 2) now provides a simple way to convert the present day value of the maximum circular 
velocity to a quantity i^peak? corresponding to the maximum circular velocity at time of 
accretion. However, we want to emphasize that this conversion has only been tested so far 
for a Milky Way sized cold dark matter halo. 

The impact of this correction scheme is illustrated in Fig. 7 for two of the simulations, 
CDM and WDM respectively. The left panels show the circular velocity curves of the 12 most 
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massive subhalos at present time {vma^{z = 0)), whereas the right panels show the sample 
when the selection is made by i^peak- Even if the selection methods do overlap to some extent, 
the right panels of Fig. 7 unsheathe quite a few subhalos that have lost a lot of their mass 
since the time of accretion. A more detailed comparison between different subhalo selection 
methods is shown in [52]. 
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